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Abstract. We performed numerical simulations of mergers between gas-rich disc galaxies, which 
result in the formation of late- type galaxies. Stars formed during the merger end up in a thick disc 
that is partially supported by velocity dispersion and has high [a/Fe] ratios at all metallicities. 
Stars formed later end up in a thin, rotationally supported disc which has lower [a/Fe] ratios. 
While the structural and kinematical properties of the merger remnants depend strongly upon 
the orbital parameters of the mergers, we find a clear chemical signature of gas-rich mergers. 

Keywords, galaxies: evolution, galaxies: formation, galaxies: interactions, galaxies: structure 



1. Introduction 

Mergers of late-type galaxies have long been associated with the formation of early- 
type galaxies (Toomre 1977). However, numerical simulations of major collisions between 
disc galaxies (e.g. Naab & Burkert 2003), which lead to the formation of ellipticals, 
have assumed a low ratio gas/stars representative of low-redshift galaxies. Going back in 
time, we expect to find fewer and fewer stars, and therefore a more important gaseous 
component. We have simulated mergers of gas-rich disc galaxies, with different mass 
ratios and orbital parameters (Brook et al. 2007; Richard et al. 2010). Our main goal is 
to find out if there is a clear observational signature of gas-rich disc-disc mergers. 

2. The Numerical Simulations 

All simulations were performed using the chemo-dynamical algorithm GCD+ (Kawata 
& Gibson 2003a, 2003b). The initial conditions consist of two galaxies with exponential 
gas discs embedded in dark matter halos. The most massive galaxy (Gall) has total mass 
of 5 X 10^ ^M© (except for simulation Mil). The mass of the second galaxy (Gal2) depends 
on the mass ratio chosen for the simulation. Fig. [T] shows the initial configuration of the 
system. Gall is at rest at Ri = (50, 0, 0) kpc. Gal2 is located at R2 = (0, 0, 15)kpc, 
with velocity V2 = (0, 100, 0)kms~^ (except for simulations M12orb and M12z). The 
baryon/DM mass fraction is 17%. Table [T] lists the parameters used for each simulation. 
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Table 1. Initials Conditions for All Simulations 



Run 




mass ratio 





V2z (kms-1) 


R2 (kpc) 


/gas 


M12 


5.0 X 10^^ 


2:1 


30.0 





0, 0, 15 


0.94 


M12orb 


5.0 X 10^^ 


2:1 


30.0 





0, 50, 15 


0.67 


M12z 


5.0 X 10^^ 


2:1 


30.0 


100 


0, 0, 15 


0.83 


M1290 


5.0 X 10^^ 


2:1 


90.0 





0, 0, 15 


0.96 


rM12 


5.0 X 10^^ 


2:1 


210.0 





0, 0, 15 


0.91 


Mil 


2.5 X 10^^ 


1:1 


30.0 





0, 0, 15 


0.88 


M13 


5.0 X 10^^ 


3:1 


30.0 





0, 0, 15 


0.92 


MHO 


5.0 X 10^^ 


10:1 


30.0 





0, 0, 15 


0.85 




Figure 1. Geometry of the initial conditions. The y-axis is pointing away from the viewer. Gall 
is initially at rest, while Gal2 is initially moving away from the viewer, in the -\-Y direction. 
The dotted lines indicate the rotation axes of the galaxies. Gall is in the XY plane, with the 
rotation axis in the Z-direction. The rotation axis of Gal2 is in the XZ plane, at an angle 
relative to the Z-axis. Both galaxies are rotating clockwise when seen from above. 

The last column lists the gas fraction at the beginning of the simulation. In each case, the 
collision leaves a very chaotic system, which then relaxes to a "quiescent state" where the 
final structure is well-established, and the star formation rate has dropped significantly. 

3. Results 

We have computed the star formation rate for all simulations (Fig. [2j). In all cases, 
a major starburst invariably occurs during the merger. After the starburst, the star 
formation steadily drops, as the system relaxes and a disc forms. We define two stellar 
populations: old stars, which are formed before and during the starburst, and young 
stars, which include all stars formed after the merger, when the starburst is completed. 

Fig. [3] show the luminosity maps of the final remnant for simulation M12. The left 
panels show young stars, the middle panels show old stars, and the right panels show all 
stars. Young and old stars form discs that have comparable radii, but the edge-on views 
(bottom panels) clearly shows that the young disc is thin, while the old disc is thick. 
Most simulations result in the formation of a thin disc made of young stars, and a thick 
disc made of old stars, except simulations M12z and Mil which show more elliptical-like 
remnants. In several cases, the merger resulted in the formation of a ring made of young 
stars. Simulations M1290 and rM12 form a central bar made of both old and young 
stars. From this analysis, we conclude that the structural properties of the remnants are 
strongly dependent of the initial conditions. 

Fig. m shows the kinematical properties for simulation M12. We calculated for each 
population the average rotational velocity V and its dispersion ay in radial bins. The 
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Figure 2. Star formation rate versus time, for simulation M12. The dashed hues indicate the 
beginning and the end of the starburst, respectively. 
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Figure 3. V-band image of M12 seen face-on (top panels) and edge-on (bottom panels), for 
young stars (left panels), old stars (middle panels), and all stars (right panels). X,Y, Z represent 
cartesian coordinates, with Z along the axis of rotation. 



left panel shows the rotation curves. In all cases, the old stars have a lower mean rotational 
velocity than the young stars. The middle panel shows {{V/cry)'^)^^'^ versus radius. In 
all cases, we find {{V/ay)'^)^^'^ > 1 for young stars, indicating the presence of a massive, 
rotationally supported disc. Most cases have ((V^/<jy)^)^/^ < 1, for old stars, indicating 
the presence of a "hot disc" supported by internal motions, but simulations M1290, rM12, 
and Ml 10 all have ((V^/<jy)^)^/^ > 1, indicating that the old stars are also rotationally 
supported. The right panel shows histograms of the stellar mass versus rotation velocity. 
In most cases, the young stars are concentrated in a narrow region of the histogram, 
while the old stars are spread in velocity. The young component is counter-rotating for 
simulation Mil, and partly counter-rotating for simulation M12z. We conclude that the 
kinematical properties of the remnants are strongly dependent of the initial conditions. 

The chemical properties are presented in Fig. [5l The top panels show [a/Fe] and 
[Fe/H] versus formation epoch for simulation M12. [a/Fe] decreases with time, while 
[Fe/H] increases strongly before and during the merger, and keeps increasing, but more 
slowly, after the merger. The bottom left panel shows [a/Fe] versus [Fe/H] for simulation 
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Figure 4. Kinematics for simulation M12. On each panel, red, blue, and black curves represent 
old stars, young stars, and all stars, respectively. Left: rotation velocity vs. radius. Middle: 
Rotational support versus radius. Right: histogram of stellar mass vs. rotation velocity, with 
negative values indicating counter-rotating stars. 




-2 -4 -2 2 



[Fe/H] [Fe/H] 

Figure 5. Top panels: a-elements abundance and metallicity versus formation epoch for simu- 
lation M12. Bottom panels: [a/Fe] versus [Fe/H] for old stars (red curve) and young stars (blue 
curve) for simulation M12 (left) and all simulations (right). 



Ml 2. Old stars have larger a-element abundance than young stars, up to a relatively 
high metallicity ([Fe/H] :^ —0.5). [a/Fe] decreases with time as the metallicity increases, 
until it reaches the values found in the thin disc. This decrease of [a/Fe] is slowed near 
the beginning of the collision, as the starburst leads to a large number of Type II SNe, 
that enrich the gas in a-elements. After the starburst. Type la SNe become effective, and 
enrich the gas in iron. This explains why [a/Fe] gradually decreases after the collision. 
The bottom right panel shows [a/Fe] versus metallicity for all 8 simulations. The results 
for old stars (red curves) are nearly identical for all simulations. Hence, the mergers 
leave a clear chemical signature: the old stars, located in the thick disc and the halo, 
have a ratio [a/Fe] that remains constant at ~ 0.5 — 0.6 with increasing [Fe/H] up to 
high metallicities ([Fe/H] = —0.5), while the young stars, located in the thin disc, have 
a lower ratio [a/Fe] which decreases with increasing metallicities. 
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